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In a quantum Markov chain, the temporal succession of states is modeled by the repeated action 
of a "bistochastic quantum operation" on the density matrix of a quantum system. Based on this 
conceptual framework, we derive some new results concerning the evolution of a quantum system, 
including its long-term behavior. Among our findings is the fact that the Cesaro limit of any 
quantum Markov chain always exists and equals the orthogonal projection of the initial state upon 
the eigenspace of the unit eigenvalue of the bistochastic quantum operation. Moreover, if the unit 
eigenvalue is the only eigenvalue on the unit circle, then the quantum Markov chain converges in 
the conventional sense to the said orthogonal projection. As a corollary, we offer a new derivation 
of the classic result describing limiting distributions of unitary quantum walks on finite graphs [J. 

PACS numbers: 03.67.-a, 03. 65 Ta, 05.40.-a 
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I. INTRODUCTION 

The theory of Markov chains, when appropriately 
generalized, provides a potent paradigm for analyzing 
the stochastic evolution of quantum systems. Over 
the past decade, motivated largely by the prospect of 
super-efficient algorithms, the theory of so-called quan- 
tum Markov chains, especially in the guise of quantum 
walks, has generated a huge volume of research, includ- 
in g m any discoveries of fundamental importance, such as 
[ll-flCj] and numerous recent advances, such as [lll - [l5j . 

In the context of quantum walks, the itinerary of the 
walker is confined to a particular topological network. 
The walker's every move, from node to adjacent node, is 
governed by a set of local rules. When applied repeatedly 
to a given initial state of the system (represented by a 
superposition of basis states) , these rules yield a succes- 
sion of new states, reflecting, ad infinitum, the evolution 
of the system. A transition rule can be either unitary 
(closed) or non-unitary (open), depending, respectively, 
on whether it is intrinsic to the system or exposes the 
system to external influences such as decoherence, noise 
or measurement. 

In this paper, we adopt the formalism of "quantum op- 
erations" , whereby both unitary and non-unitary rules of 
state transition, as well as various combinations thereof, 
are treated under a unified mathematical model. In 
this framework, the "transition matrix" of a classical 
Markov chain is replaced by a "bistochastic quantum op- 
eration" and the "state distribution vector" of the classi- 
cal Markov chain is replaced by a "density matrix" . The 
resulting description of quantum state evolution, known 
as a quantum Markov chain [J, [HI , turns out to resemble 
very closely the evolution of a classical Markov chain. In 
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particular, the same kind of questions pertain, including 
questions about the long-term evolution of the process 
and the possible existence of a limiting state. 

Among our findings is the fact that the Cesaro limit 
of any quantum Markov chain converges always to a sta- 
tionary "state" , regardless of the initial state. As a note- 
worthy special case of this result, we remark that for any 
unitary quantum walk on a graph, as in the limit of 
the time-averaged probability distribution always exists. 

To complete the picture, we specify conditions for the 
existence of a limiting state in the strict (non-Cesaro) 
sense of the word "limiting" . In the strict sense, it turns 
out that the limiting behavior depends only on the de- 
ployment on the unit disc of the eigenvalues of the bis- 
tochastic quantum operation. Specifically, if A = 1 is the 
only eigenvalue on the unit circle, then, for any given 
initial state, the associated quantum Markov chain con- 
verges to a stationary state. Moreover, if the eigenspace 
of A = 1 is one-dimensional or contains only a single 
density operator, then the associated Markov chain con- 
verges to the maximally mixed state, irrespective of the 
initial state. Otherwise, if the bistochastic quantum op- 
eration possesses any eigenvalue on the unit circle other 
than A = 1, then a limiting state is not guaranteed to 
exist except in the generalized sense of Cesaro. These 
findings are seen to be analogous, in a very natural way, 
to the fundamental properties of classical Markov chains 
(e.g. in [17|, Chapter 8). 

Our results may represent substantial progress toward 
answering the first of a set of "open questions" posed 
by Ambanis in 2005 [lj|. We would be remiss not to 
acknowledge our indebtedness to which treats the 
special case of a quantum Markov chain generated by a 
random unitary operation. 

In what follows, after a brief review of some prelimi- 
naries (Section 2), we proceed (Section 3) to present our 
main results on the question of the limiting behavior of a 
quantum Markov chain. The proofs of the results given 
in Section 3 are deferred to Section 4. In Section 5, we 
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present a comprehensive classification of quantum oper- 
ations according to their limiting behavior. Finally, in 
Section 6, we offer some concluding remarks, including 
some relevant questions for further investigation. 



II. BISTOCHASTIC QUANTUM OPERATIONS 

Given a Hilbert space H of finite dimension N, let 
03(H) denote the set of all linear operators on H, with 
inner product defined by 

(X,Y)=tr(XlY). (1) 

The corresponding norm, called Frobenius norm or 
Schatten 2- norm, is defined by 

||X|| = [trpftx)] 1 ^ = v^X). (2) 

This choice of norm on 05(H) will remain in effect 
throughout this paper. 

Let D (H) C 05(H) denote the set of positive operators 
p-.U^rU with Tr(p) = 1. The operators p G D(H) are 
the so-called "density operators". They serve to model, 
as faithfully as do the "state vectors" themselves, the 
possible states of a quantum system whose state vectors 
reside in H. 

By a super- operator <fr on 05(H), we mean a linear 
mapping : 55(H) — > 55(H), with norm defined by 

11*11 ^Su Pxe2HW) M^. (3) 

Note that dim05(H) = N 2 , where N = dim(H). Thus, 
any super-operator on 05(H) can be represented, rel- 
ative to a given basis for 05(H), by an TV 2 x iV 2 matrix. 
In the sequel, this matrix will be denoted by the sym- 
bol [<&]. In particular, relative to a special basis con- 
sisting of eigenvectors and generalized eigenvectors of 
the shape of the matrix [<&] conforms to a special quasi- 
diagonal lay-out called the Jordan canonical form. The 
details can be found in any one of a number of sources, 
including [l8j . 

Among the set of super-operators, we distinguish a spe- 
cial subset called "quantum operations". By definition, 
to qualify as a quantum operation, the super-operator $ 
must be completely positive, meaning that the extended 
map 3» ® I„ is positive for all n > 1 . 

The formalism of quantum operations is flexible 
enough to handle both unitary (closed) and non-unitary 
(open), or a mixture thereof, of discrete transitions of 
state of a quantum system. For a g ood introductory ex- 
position of this subject, see 0, |20|| 

By Choi's Theorem [U and [ll, SS SI, any com- 
pletely positive linear operator, including any quantum 
operation <& : 05(H) —> 05(H), can be represented in 
terms of a set A = {Ai \ i = 1, 2, iV 2 } of "Kraus oper- 
ators" , as follows: 



* A {X)=Y,AiXA\. (4) 

i 

In this expression, which we call the "Choi expansion" of 
€>, the symbol A\ denotes A^ (transpose of the complex 
conjugate of Ai). 

In terms of the Choi expansion, the condition of being 
trace-preserving, meaning that Tr(3?_4(A)) = Tr(A) for 
all X G 05(H), is equivalent to the condition: 

E4A = %. (5) 

i 

On the other hand, if the Kraus operators of satisfy 
the dual condition: 

£>4=%, (6) 

i 

then <1?.4 is said to be unital. Note that (6) is equivalent 
to the simple statement that $^(Ia?) = Ijv- 

A quantum operation which is both unital and trace- 
preserving is called bistochastic. The term doubly 
stochastic quantum channel also has been used to refer 
to quantum operations of this sort [231 ]. Note that a 
bistochastic quantum operation transforms elements of 
£>(H) into elements of S(H). In other words, since the 
elements of 23(H) represent the states of a quantum sys- 
tem, a bistochastic quantum operation transforms states 
of a quantum system into other legitimate states of that 
system. 

III. LIMIT THEOREMS FOR BISTOCHASTIC 
QUANTUM OPERATIONS 

In the sequel, the proofs of all theorems, corollaries 
and supporting lemmas are deferred to Section 4. 

By [24[, a bistochastic quantum operation <fr_4 must 
satisfy the condition \\&a\\ < 1- Thus the spectrum of 
<fr_4 is confined to the unit disk. Also, since *^(Ijv) = 
In, we see that A = 1 is an eigenvalue of 3>_4 and \\&a\\ — 
1. For future reference, we record these observations in 
the form of a lemma: 

Lemma 1. Let <I>_4 be a bistochastic quantum opera- 
tion on the Hilbert space 05(H), then 

i- II*aII = i 

2. If A is an eigenvalue of <fr.4, then |A| < 1. 

3. The value A — 1 is an eigenvalue of 3? ,4. 

The observations recorded in Lemma 1 are new by no 
means. See, for instance, [25| . 
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For an eigenvalue A of let Ker(_*_4 — AI) and 

Ran(<fr_4 — AI) denote, respectively, the kernel and range 
of the operator — AI on the Hilbert space 58(H). 

Lemma 2. Let be a bistochastic quantum op- 
eration on 58 (H) and let A be an eigenvalue of &a with 
|A| = 1. Then Ker(*^ - AI) n Ran(*^ - AI) = {0}. 

From the preceding lemma, we can derive an important 
inference concerning the algebraic and geometric multi- 
plicities of the eigenvalues of $>a- For an eigenvalue A 
of <J?_4, let its algebraic multiplicity be denoted by m(A) 
and lets its geometric multiplicity be denoted by 17(A). 
Recall that g(X) = dimKer( _>_4 — AI). In addition, let the 
spectrum of 3? .4 be denoted by A(<__4) and let Ai($_^) 
denote the subset of A( <_^) consisting of A G A( <__4) with 
|A| = 1. 

Lemma 3. If A G Ai( _>_4), where is a bistochastic 
quantum operation on 03(H), then m(A) = g(A). 

As in the proof of Lemma 3, let [&a] denote the 
TV 2 x A^ 2 Jordan canonical matrix representation of <&a- 
To conserve type-set space, and without undue risk of 
confusion to the reader, we prefer to display the matrix 
[_\a] m the following self-explanatory format: 

[&a] = diag(Ai, A 2 , A fc , J u J 2 ,-, A), (7) 

where Aj G Ai(<__4) and J r , r = 1,2,..., ft,, denote the 
Jordan blocks corresponding to eigenvalues whose norms 
are strictly less than unity. 

For an eigenvalue A of 3>_4, let E& A (A) = Ker( <__4 — AI) 
denote the eigenspace of A. In particular, for A = 1, it 
has been established [26-28] that 

£ # ^(l) = Ker(*4-I) 
= {X G 05(H) : XA t = A t X;i = 1,2,...}. (8) 

The following lemma articulates the special status of 
E& A (1) relative to the other eigenspaces of <_^. 

Lemma 4. Let A be an eigenvalue of €>_4 such that 
I A I = 1 and A 7^ 1. Let a be an eigenvalue of with 
a I < 1. Let Yi, Yj a denote the generalized eigenvec- 
tors belonging to a. Then 

1. ^(l)l^(A) 

2. E* A (1) ±Spm{Y u ...,Y ja } 

In other words, by Lemma 4, the eigenspace E<f, A (l) is 
orthogonal to every eigenvector (including every gener- 
alized eigenvector) of belonging to eigenvalues other 
than A = 1. 

Theorem 5. Let be a bistochastic quantum op- 
eration on the Hilbert space 03(H) and let p(0) G 55(H) 
denote the density matrix representing the initial state of 
a quantum system. Let g(l) denote the geometric mul- 
tiplicity of the eigenvalue A = 1 and let {^r}i<r<g(i) 



denote an orthonormal basis for E^ A (1). If A = 1 is the 
only eigenvalue of &a on the unit circle, then the iterated 
succession of quantum states p(t) = & A p(0) converges to 
p(oo) = J2 9 i=^ T ( z jp(0))Zi. In particular, if g(l) = 1, 
then lim.t_j.oo $> A p(0) = -^1, independently of the initial 
state p(0). 

In previous publications, such as [Hj], the special case 
of a bistochastic quantum operation whereby g(l) = 1 
is known as a primitive quantum channel. The reference 
[29l | provides an excellent analysis of this line of investi- 
gation. 

In the literature, numerous examples can be found 
of bistochastic quantum operations (a.k.a. quantum 
channels) for which the limiting behavior of the asso- 
ciated quantum Markov chain is governed by Theorem 
5. Among these examples are 1) the quantum walk on a 
cycle, subject to decoherence on the degree freedom of 
coin [30| and 2) the generalized depolarizing channel on 
a multiple-dimensional quantum system [lj| [2(| . These 
examples can be treated as special cases of the following 
corollary to Theorem 5. 

Corollary 6. If the bistochastic quantum operation 
&A is given by 

^ A {X) = {l-p)X + Y J A i XA\, (9) 

i 

where ^z^j = pi and < p < 1, then, as in Theo- 
rem 5, p(oo) = £r=i tr(Z}p(0))Z r . 

In the context of quantum channels, the Kraus opera- 
tors defining the generalized depolarizing channel can be 
scalar multiples of the discrete Weyl operators (or the 
generalized Pauli operators). When such is the case, as 
in for instance [23j, the eigenspace of the eigenvalue 1 
must be one-dimensional. Consequently the associated 
quantum Markov chain must converge to the maximally 
mixed state, regardless of the initial state. 

The main assertion of Theorem 5 pertains only to bis- 
tochastic quantum operations possessing no eigen- 
values on the unit circle other than A = 1. In general, 
when this condition is relaxed to admit other eigenvalues 
on the unit circle, the expression p(t) = $ A p(0) is com- 
pelled no longer to converge to a stationary state. Many 
examples can be found in the literature, including classic 
cases of unitary quantum walks on finite graphs [l| . 

However, even if a limiting state fails to exist in the 
usual sense, we still might want to probe the possibility 
of a "limiting state" p(oo) in the sense of Cesaro: 

= lim ^(o) + ^(o) + - + ^p(o), (10) 

t-voo t 

In terms of this generalized sense of "limiting state" , it 
turns out that every quantum Markov chain converges. 
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Theorem 7. Let be a bistochastic quantum op- 
eration on the Hilbert space 25 (H) and let p(0) G D(H) 
denote the density matrix representing the initial state 
of a quantum system. Let g(l) denote the geometric 
multiplicity of the eigenvalue A = 1 and let {Z r } 1<r<g ^ 
denote an orthonormal basis for Eq, A (1). Then, for every 
p(O)g0(H): 



p(oo) 



9(1) 

J2tr(zjp(0))Zi. (11) 
i=i 



According to this theorem, the sequence of Cesaro 
means {j Y^n=i ^j\}t=i ^ s guaranteed always to converge 
to a limit which coincides with the orthogonal projec- 
tion of the initial state upon the eigenspace of the unit 
eigenvalue of the bistochastic quantum operation. It is 
only fair to point out that alternate versions of this re- 
sult have appeared in the literature. For instance, in 
[26j | (see Theorem 2.4(a)) the projection on the invariant 
subspace is obtained as a weak limit of a subsequence 
of the full sequence of Cesaro means. When applied to 
the initial state, this projection yields the desired fixed 
point. Another close match to Theorem 7, derived by 
different means and framed in different language, can be 
found in Chapter 6 of [3l|. However, to the best of our 
knowledge, none of the published variants of this theo- 
rem have explicitly inferred that the projection must be 
orthogonal. 

As an immediate corollary of Theorem 7, it follows, 
as in [l| , that the time-averaged probability distribution 
for any unitary quantum walk on a finite graph must 
converge. For a unitary quantum walk on a finite graph 
starting with an initial state |ao), let P n (v\ao) denote 
the probability of finding the walker at node v at time 
n, and let P t (v\ao) = | J2n=i Pn(v\ao) denote the time- 
averaged probability. 



Corollary 8. (see Theorem 3.4., [l|) Let ipj, Xj de- 
note the unit eigenvectors and corresponding eigenval- 
ues of the unitary operator U associated with a unitary 
quantum walk on a finite graph and let the bistochastic 
quantum operation <1>_4 be defined by ^a(p) = UpW. 
Then, for any initial state |ao) = J2j a j\^j)i we have 

lim P t (v\a ) = V* a l a*(a,v\tp i )(iljj\a,v}, (12) 

£— s-oo * — * J 
i,j,a 

where the sum is restricted to pairs i, j such that A^ = Xj. 

IV. PROOFS OF THE THEOREMS 

This section is reserved for the proofs of theorems, lem- 
mas and corollaries given in Section 3. 



A. Proof of Lemma 2 

Suppose X G Ker(<f> A - XI) n Ran($^ - AI). Then 
<f> A (X) = XX and &a(Y) - XY = X for some Y G 
55(H). Applying the linearity of &a, we infer that 
(Y) = X n Y + nX^X for all n > 1. Conse- 
quently WnX^XW = \\$ A (Y) - A"F||, which implies 
that n\\X\\ < + But since = 1 (see 

Lemma 1), we have n||X|| < 2||Y||. Since this inequal- 
ity must hold for all n > 1, we conclude that \\X\\ = 0, 
whence X = 0. □ 

In the above proof, we have borrowed liberally from 
the reasoning employed in [l8j], which treats the special 
case of a random unitary operation. 



B. Proof of Lemma 3 

We proceed by contradiction. Suppose m(X) ^ 5(A). 
Let [^a] denote the N 2 x N 2 matrix representation of 
&A i n Jordan canonical form. On the one hand, [&a] 
must contain a Jordan block J belonging to A of size 
> 1. On the other hand, by standard matrix theory, 
there must exist a generalized eigenvector, say v, of &a 
such that (4>_4 — XT)v is itself an eigenvector of 3>a- This 
implies that Ker(& A - AI) n Ran(& A - AI) ^ {0}, which 
contradicts Lemma 2. □ 



C. Proof of Lemma 4 

The adjoint operator of <&a is given by 

& A {X) = Y,A\XA l . (13) 



Thus, by Eq. ©, we have Ker(*^ - I) = Ker(*^ - I). 

To prove statement 1, let Z G Ker(€>_4 — I) and Y G 
Ker(*4 - AI). Then (Z,Y) = ($ A Z,Y) = {Z,& A Y) = 
X(Z,Y). Since A ^ 1, it follows that (Z,Y) = 0. Hence 
Ker(*A - I) _L Ker(*4 - AI). 

We proceed to justify statement 2. For an eigenvalue 
a with |a| < 1, we may assume, without loss of general- 
ity, that the generalized eigenvectors belonging to a are 
arranged in a sequence Yi, Y%, Yj a such that: 



{&A ~ al)Y r = y r _!, 



(14) 



where, by definition, Yq = 0. It follows that (&a — 
al) r Y r = 0. If Z G Ker(*^ - I), then (Z,Y t ) = 
{^ A Z,Y X ) = (Z,$aYi) = a{Z,Y t ), which implies 
that (Z,Yi) = 0. Similarly, (Z,Y 2 ) = (<f> A Z,Y 2 ) = 
(Z,$ A Y 2 ) = (Z,Yi) + a{Z,Y 2 ), from which it follows 
that (Z, Y 2 ) — 0. Likewise, by the same reasoning, ap- 
plied repeatedly, we deduce that (Z,Y r ) = for any r. 
Thus Ker(*A-I) _L Span{Yi, Y 2 , Y ja }. □ 
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D. Proof of Theorem 5 



Proof of Theorem 7 



Let «i, ah denote the eigenvalues of &a whose abso- 
lute values are less than 1. Since, by Lemma 3, A = 1 is 
the only eigenvalue on the unit circle, a basis for 03 (Ti.) 
might be assembled as follows. Starting with an or- 
thonormal basis for eigenspace E^ A (1), append a maxi- 
mal set of linearly independent generalized eigenvectors 
belonging to the eigenvalues a ri r = 1,2,..., ft,. Then, 
relative to this basis, the Jordan canonical matrix repre- 
sentation of is given by: 



diag I 



Jh) 



(15) 



where J r is the Jordan block corresponding to the eigen- 
value a r , r = 1, 2, h. 

Consider what becomes of the Jordan blocks of the 
powers [^a]* as t — > oo. Since each of the Jordan blocks 
J r is an upper triangular matrix whose diagonal is popu- 
lated by a single eigenvalue of modulus strictly less than 
unity, it is a simple exercise in elementary algebra to show 
that lim.t_j.oo J* = O r (zero matrix of same size as J r ). 
Thus, if we define 



=diag (I fl(1) ,0 1 ,0 3 ,...,0 h ) 



(16) 



then \\3>\ - $3" || -> 0. 

Next, we consider the effect upon on the initial state 
p(0) of & A as t — > oo. 

By Lemma 4, we may choose for 03(H) a basis con- 
sisting of a basis Z\, Z2, ..Z g i\\ for E& A (1) together with 
an orthogonally complement basis consisting of general- 
ized eigenvectors belonging to all other eigenvalues. In 

W 



terms of such a basis we have p(0) — E?=i c rZ r 
where W _L Z r for all r = 1, 2, ...,g(l). By sim- 
ple linear algebra, it follows that c r = tr(_?Tp(Q)) for 
r = 1, 2, ...,g(l). Therefore, by Eq. (fll)|) . we have 

lim^oo *>(0) = J2 9 Sl tr{Z$ P {Q))Z r . 

Finally, if dimKer( — I) = 1, then, to serve as a nor- 
malized basis for E^ A (1), we may take {-7=!^}, from 

which it follows that limt_>. 0O <& A p(0) = -^-1, for any ini- 
tial state p(0). □ 



To evaluate the limit in Eq. (fTO|) . we have only to re- 
examine Eq.([7]). For every eigenvalue A = A^ of <__4, 

>t m 1 A-A t+1 



A" 



if A ^ 1, then |£n=l" - nix- 
larly, for every Jordan block J = J r in Eq.©, we have 
7EL1 J " = K 1 _ J )~V ~ ^* +1 ) ~> 0. Based on 
these two observations, we deduce that 7 
converges to the diagonal matrix 



— ► 0. Simi- 



Ml) 











(17) 



The expression for p(oo) follows by a path of reasoning 
identical to that employed in the proof of Theorem 5 to 
obtain the analogous expression for the limiting state in 
that situation. □ 



G. Proof of Corollary 8 

As an orthonormal basis for the eigenspace E& A (1), wc 
may take the set of vectors \J i jilipi/i^j |}> where the 
union is restricted to pairs i,j such that Xi — Xj. Ac- 
cordingly, the orthogonal projection of the initial state 
p(0) = \ao)(ao\ into the eigenspace E^, A (1) is 

£(|^<V_l,p(0))hM(^l =£a* a ihM<^l. (18) 



where, as above, the sum is restricted to pairs i,j such 
at Xi - 
Since 



that Xi = Xj 



Pt(v\a ) = ^tr(\a,v)(a,v\^^\a )(a \\ , (19) 



and since tr (|a, v )(a, v\-) is a continuous function of 
the argument, it follows, by Theorem 7, that Pt(v\ao) 
converges to 



E. Proof of Corollary 6 



J^tr \a,v){a,v\^a l a*\ip l ){ip j \ =^aia* j (a,v\i> i )ty j \a,v)2ti) 



• ■j 



It suffices to verify that A = 1 is the only eigenvalue of 
$>A on the unit circle. By (24|, whenever ^ AiA\ = pi, 

then HEi^^ill < P\\ X \\ for an Y operator X. Now 
suppose A is an eigenvalue of on the unit circle and 
X is an eigenvector belonging to A. Then _»_4(X) = 
XX = (1 - p)X + J2i A t XA\, from which it follows that 
||AX|| = \\X\\ = \\(l-p)X + j: i A i XAl\\ < ||(l-p)__|| + 

||Ei^4ll.< (1 - P)\\X\\ + P\\*\\ = 11*11- Thus, all 
of the preceding inequalities actually are equalities. In 



particular, we have || J2i AiXAl\\ = p\\X\ 
p)X — J2i AiXA], which implies that |A 
p\\X\\. It follows that A = 1. □ 



Thus (A-l- 
l+p|-||A|| = 



where the sum is restricted to pairs i, j such that A^ = Xj. 
□ 



V. CLASSIFICATION OF BISTOCHASTIC 
QUANTUM OPERATIONS 

Let &a ■ S(H) — > 23(H) be a bistochastic quantum 
operation and let p(0) _ 55(H) denote the density ma- 
trix representing the initial state of a quantum system. 
Depending on the limiting behavior, as t 00, of the cor- 
responding quantum Markov process, <&a must belong to 
one of the following four mutually exclusive categories: 
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(1) lim t _ i . 00 <fc A p(0) converges to the maximally-mixed 
state jjl, independently of the initial state p(0). 

(2) lirrii-j.oo & A p(0) converges, but the limit depends 
upon the initial state p(0). 

(3) lim^oo & A p(0) fails to converge, but the Cesaro 
average limt_j.oo \ J2 n =i ^aP(®) exists an d equals the 
maximally-mixed state -^1, independently of the initial 
state p(0). 

(4) Hindoo <f> f A p(0) fails to converge, but the Cesaro 
average lim.t_j.oo | Yl n =i ^aP(®) exists and depends upon 
the initial state p(0). 

To elucidate each of the above categories, we proceed 
to offer some comments and examples. 

In category (1), the quantum operation _»_4 possesses 
no eigenvalue on the unit circle other than A = 1. 
Moreover, the eigenspace Ker( _>_4 — I) of A = 1 is one- 
dimensional, spanned only by the density operator -^1, 
or contains only a single density operator -^1. 

As an example, consider the quantum operation [l6| 
defined by 



* A (X) = -X 



1 



uxtf 



where U is the unitary transformation given by 



U 



1 __ 

2 

V3 1 

2 2 



(21) 



(22) 



According to 18], this is an example of a so-called 
random unitary operation. The eigenspace Ker( _»^ — AI) 
where A _ Ai(<_^) is equal to the set D\ := {X e *B(H) : 
IX = XXI and UX = XXU}. Evidently, Ai(*a) = {1}. 
A simple calculation shows that the eigenspace Ker( &a 
I) of the eigenvalue A = 1 is D\ = {kl : k 6 C}, which is 
one-dimensional. Therefore lim.t_j.oo & A X — 5I2 for any 
initial state X. 

It can be verified that the bit-phase flip channel (p377 
in [19() also belongs to category (1). A less trivial exam- 
ple of this sort can be found in [32j |. 

In category (2), possesses as its only eigenvalue on 
the unit circle the value A = 1 and Ker( — I) is at least 
two-dimensional. 

As an example, consider the quantum operation _>_4 
associated with a so-called phase flip channel on single 
qubits [li|, given by 



$ A (X) =pX + (l- p)ZXZ^ 



where Z is the Pauli matrix: 



1 

-1 



(23) 



(24) 



By a pattern of reasoning similar to that employed in 
the previous example, we infer that Ai(<__4) = {1}, and 
the eigenspace Ker( — I) is spanned by the two linearly 
independent density operators: 



X 1 



1 


0" 


, x 2 = 


"0 


0" 











1 



(25) 



In this case, lim t _j. 
are given by: 



00 ™ A 



Xoo where X and A 



X = 



an 0,12 

0-21 0-22 



Aoo — 



an 
a 2 2 



(26) 



It can be verified without much difficulty that the bit 
flip channel ([Hj], P-376) belongs also to category (2). 
Additional examples of this sort can be found in ([_0j. 
p.108). 

In category (3), the quantum operation & A possesses 
at least two distinct eigenvalues (including A = 1) on the 
unit circle. The eigenspace of A = 1, namely Ker(_>^— I), 
is one-dimensional and spanned by I. An example of this 
type of quantum operation is provided by quantum walks 
on the TV-cycle. In this scenario, N is assumed even, the 
Hadamard transform serves as the coin operator and the 
evolution of the system is subject to decoherence on both 
the position and the coin degrees freedom. A detailed 
treatment of this model is planned for our forthcoming 
paper. 

In category (4), has at least two distinct eigenval- 
ues (including A = 1) on the unit circle and Ker( _>_4 — I) 
contains at least two linearly independent density 
operators. As an example of this type of quantum 
operation, we cite [18[ which studies the properties 
so-called two-qubit controlled-not operators. 



VI. CONCLUSION AND RELATED QUESTIONS 

Evidently, for a stochastic quantum operation <__4, the 
eigenvalues lying on the unit circle determine the evo- 
lution of the associated quantum Markov process, in- 
cluding the existence or non-existence of a long-term sta- 
tionary state. More precisely, the long-term behavior of 
the quantum Markov process is intimately linked to the 
structure of the eigenspaces of eigenvalues on the unit 
circle. 

We speculate that the eigenspace of an eigenvalue A on 
the unit circle, denoted by D\, might conform always to 
a formulation in terms of a set A = {Ai | i = 1, 2, N 2 } 
of "Kraus operators" , as follows: 

D\ = {X G SCH) : XAi = XAiXiori = 1,2,...}. (27) 

If valid, this formulation conceivably might provide an ef- 
ficient means for identifying the eigenspaces of all eigen- 
values of absolute value 1. The origin of this idea can be 
attributed to [H, [H-ill . 
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